Ni-substituted sites and the effect on Cu electron spin dynamics of YBa 2 Cu 3 - x Ni x 7 _s 
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We report Cu nuclear quadrupole resonance experiment on 
magnetic impurity Ni-substituted YBa2Cu3_ :E Ni :E 07-_,5. The 
distribution of Ni-substituted sites and its effect on the Cu 
electron spin dynamics are investigated. Two samples with 
the same Ni concentration x=0.10 and nearly the same oxy- 
gen content but different T c 's were prepared: One is an as- 
synthesized sample (7-5=6.93) in air (T c ~ 80K), and the 
other is a quenched one (7-<5=6.92) in a reduced oxygen at- 
mosphere (T c ~ 70K). The plane-site 63 Cu(2) nuclear spin- 
lattice relaxation for the quenched sample was faster than 
that for the as- synthesized sample, in contrast to the 63 Cu(l) 
relaxation that was faster for the as-synthesized sample. This 
indicates that the density of plane-site Ni(2) is higher in the 
quenched sample, contrary to the chain-site Ni(l) density 
which is lower in the quenched sample. From the analysis 
in terms of the Ni-induced nuclear spin-lattice relaxation, we 
suggest that the primary origin of suppression of T c is as- 
sociated with nonmagnetic depairing effect of the plane-site 
Ni(2). 
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Magnetic impurity causes depairing in both s-wave and 
d-wave superconductivity. Complete suppression of the 
superconducting transition temperature T c is observed 
for Ni-doped I^-xSr^CuO^ and YBa2Cu40s, somewhat 
more weakly than for Zn-doped ones Ni impu- 

rity in the high-T c cuprate superconductor carries a lo- 
calized moment, because the uniform spin susceptibility 
with Curie or Curie- Weiss law is observed. The depair- 
ing effect of potential scatterer Zn on the d-wave super- 
conductivity is a natural consequence from breakdown 
of Anderson's theorem For YBa2Cu 3 07_(5 with 

the optimized T c =92 K, however, the decrease of T c per 
Ni concentration is smaller than that per Zn concentra- 
tion ||-|| . Figure 1 (a) shows the impurity doping depen- 
dences of T c for various high-T c superconductors with 
Ni or Zn The decrease of T c by Ni doping for 

YBa2Cu307_,5 is smallest among these compounds with 
Ni in Fig. 1(a). This small decrease of T c by Ni doping 
for YBa2Cu307_,5 was attributed to a weak scattering of 
conducting carriers by Ni impurities (Born scatterer) |l3| ] 
or to a softening of pairing frequency itself |9| . Neverthe- 
less, it has been suspected that only a part of the doped 
Ni impurities is substituted for the plane Cu(2) site and 
that the remaining part is for the chain Cu(l) site 14 1?]] 



(see the references in Ref . [Q ) . The bulk T c is considered 
to be determined by the amount of the in-plane Ni(2) 
impurity. One should note that YBa2Cu307_i has two 
crystallographic Cu sites, i.e. the chain Cu(l) and the 
plane Cu(2) sites. The recent observation of an in-plane 
anisotropy of optical conductivity for detwinned single 
crystal YBa2Cu3_ a; Ni 2 ;07_ ( 5 indicates the existence of the 
chain-site Ni(l) p9| . A microscopic study on the selec- 
tive substitution of Ni impurity is therefore of interest. 

Synthesis of oxides under reduced oxygen partial 
pressure is frequently effective in controlling cation 
solid solution, e.g. to synthesize the superconduct- 
ing Lai +a; Ba2-a:Cu307_5 p0|] , or to optimize the su- 
perconducting critical current density or the irreversible 
magnetic field of Ndi +a; Ba2_ x Cu307_5 21 2g]. Re- 
cently, using the reduced oxygen partial pressure tech- 
nique, Adachi et al., succeeded in controlling T c of 
YBa2Cu3_ 2; Ni 2 ;07_5 with optimal oxygen content [ pj| . 
They synthesized two series of YBa2Cu3_ 2; Ni 2 ;07_5 sam- 
ples having the different T c 's per Ni concentration. Fig- 
ure 1(b) shows the Ni doping dependence of T c of their 
samples Q. Here, we focus on two samples with the 
same Ni concentration x—0.10 (2=0.033 in Fig. 1(b)) 
and nearly the same oxygen content but different T c . 
For convenience, let us call one sample with T c « 80K 
(7-5=6.93) an as-synthesized one, because it was syn- 
thesized in flowing oxygen gas without quenching treat- 
ment, and the other sample with T c « 70K (7-5=6.92) a 
quenched one, because it was the as-synthesized sample 
once again fired and quenched in a reduced oxygen atmo- 
sphere at 800 °C. The details are given in Ref. pi Ni 
prefers the higher coordination of oxygen atoms~p3|,f7| . 
The plane-site Cu(2) is located in the pyramid with five 
oxygen ions, whereas the chain-site Cu(l) is coordinated 
with two, three, or four nearest neighbor oxygen ions. 
The two series of samples with different T c suggest that 
the distribution of Ni-substituted sites over Cu(l) and 
Cu(2) sites is changed through synthesis under the re- 
duced oxygen atmosphere. 

In this paper, we report the measurements of Cu(l) 
and Cu(2) nuclear quadrupole resonance (NQR) spec- 
tra and nuclear spin-lattice relaxation curves, to study 
microscopically the distribution of Ni-substituted sites 
for the above mentioned two samples (as-synthesized and 
quenched ones). The observed difference in the 63 Cu nu- 
clear spin-lattice relaxation at Cu(l) and Cu(2) in the 
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two samples must be a Ni doping effect. From the Cu 
NQR measurements, we give proof that the heat treat- 
ment in reduced oxygen atmosphere results in a redistri- 
bution of the Ni atoms in YBa2Cu3_ x Ni I 07_i. 

Two samples with precisely the same mole number 
could not be prepared, because some parts of the sam- 
ples have already been used for the characterization [El| . 
Hence, quantitative comparison of the relative intensity 
of Cu NQR spectra could not be made to estimate the 
relative number of the observed nuclei. The nuclear spin- 
lattice relaxation is independent of the sample volume 
and is relatively more sensitive to the impurity than the 
intensity of NQR spectrum. The powder samples were 
coated in paraffin oil. A coherent- type pulsed spectrome- 
ter was utilized for the zero-field Cu NQR measurements. 
The Cu NQR frequency spectra with quadrature detec- 
tion were measured by integration of the spin-echoes as a 
function of rf frequency. The Cu nuclear spin-lattice re- 
laxation curves were measured by an inversion recovery 
spin-echo method, where the Cu nuclear spin-echo am- 
plitude M{t) was recorded as a function of time interval 
t after an inversion 7r pulse, in a 7r — t — n/2 — 7r-echo 
sequence, and M(oo) was also recorded in a tt/2 — 7r-echo 
sequence (no inversion pulse) as usual |ll],[l^,|4|,|5| . 

Figure 2 shows the Cu NQR spectra at T=4.2 K. The 
amount of the as-synthesized sample is more than that 
of the quenched one, nevertheless the intensity of Cu(l) 
NQR spectra for the as-synthesized sample was weaker 
than that for the quenched one. But, we could not find 
qualitative difference in the line profiles of Cu NQR spec- 
tra between two samples. 

Figure 3 shows the 63 Cu nuclear spin-echo recov- 
ery curves (spin-lattice relaxation curves) pit) = 1 — 
M(t)/M(oo) for Ni-doped samples at T=4.2 K. Inset fig- 
ures show the recovery curves for pure (impurity-free) 
YBa2Cu30e.98 (T c =92 K). Solid curves are the least- 
squares fitting results using the theoretical function de- 
scribed below. First, from comparison with inset fig- 
ures, all the recovery curves for both Ni-doped samples 
recover more quickly than those for pure YBa2Cu3 06.98- 
Thus, Ni impurities distribute over both sites of Cu(l) 
and Cu(2). Second, the Cu(2) nuclear spin-echo signal 
recovers faster in the quenched sample than in the as- 
synthesized one, whereas the Cu(l) nuclear spin-echo sig- 
nal recovers slower in the former than in the latter. The 
Cu(2) nuclear spin-echo signal is affected in the quenched 
sample more than in the as-synthesized one, whereas the 
Cu(l) is vice versa. Thus, it is natural to conclude that 
the amount of Ni(2) in the quenched sample is more than 
that in the as-synthesized one and that of N(l) is vice 
versa. 

The 63 Cu nuclear spin-echo recovery curves for Ni- 
doped samples in Fig. 3 are nonexponential func- 
tions. For quantitative discussion, we analyzed the 
experimental recovery curve p{t) by the exponential 
function times a stretched exponential function p(t) = 



p(0)exp[- w t/(T 1 ) HO 5T - yjwt/n] (p(0), {Tx)host and 
Ti are the fit parameters, and w—3 at Cu(2) and 
w=l at Cu(l)), after the dilute magnetic alloy |^6| or 
YBa2Cu408 with Ni impurities |24|]. The multiplicative 
numerical factor w is introduced to conform to the con- 
ventional expression of T\ [ p7| , and it is not essential in 
the below discussion. Here, w — 3 is defined for the 
Cu(2) NQR T\ under a uniaxial electric field gradient. 
The Cu(l) site is under an asymmetric electric field gra- 
dient p8| , p9] |, so that all the x-, y-, and z-components of 
the fluctuating local field contribute to the Cu(l) NQR 
T\. Then, we use a simple w=l. (Ti)host is the Cu 
nuclear spin-lattice relaxation time due to the host Cu 
electron spin fluctuation via a hyperfine coupling. T\ is 
the impurity-induced nuclear spin-lattice relaxation time 
via a longitudinal direct dipole coupling or a two dimen- 
sional Ruderman-Kittel-Kasuya-Yosida interaction f{0|| . 
We have confirmed no significant contribution from nu- 
clear spin diffusion by measuring 65 Cu isotope depen- 
dence and pulse-strength H\ dependence of the recovery 
curves |31 . Although for the heavily impurity-doped sys- 
tem |}2 1 or spin glass system j3j| where it is hard to assign 
separately the host and the guest contributions, a single 
stretched exponential function p(t) = p(0)exp[— (i/ T l) a ] 
with a variable exponent a might be appropriate, we be- 
lieve that the present model with two time constants, 
(Ti)host and n, is minimal and appropriate for the 
x=0.10 samples. 

Figure 4 shows the temperature dependence of the es- 
timated {l/Tx) HO ST (a) [(c)] and of 1/n (b) [(d)] at 
63 Cu(2) [ 63 Cu(l)]. Below T c , the Ni-induced relaxation 
component (stretched exponential part) is predominant 
both at Cu(l) and Cu(2). In Figs. 4(b) and 4(d), the Ni- 
induced relaxation rate 1/n of Cu(2) for the quenched 
sample is more enhanced than that for the as-synthesized 
one, whereas that of Cu(l) is vice versa. In Fig. 4(b), the 
upturn of 1/n of Cu(2) below 10 K for the as-synthesized 
sample is consistent with the upturn of an initial relax- 
ation rate 1/T is reported in Ref. jjJJ. In Fig- 4(d), the 
difference in 1/ti of Cu(l) between the two samples with 
Ni is larger at lower temperatures than about 30 K. It is 
hard to estimate precisely the small 1/ti and its differ- 
ence above about 30 K, where the signal-to-noise ratio of 
Cu(l) NQR is poor in the Ni-doped samples. 

In Fig. 4(a), the host relaxation rate (1/Ti)host of 
Cu(2) in the as-synthesized and the quenched samples 
decreases more steeply than that in pure YBa2Cu30g.98, 
as the temperature is decreased below T= 18^20 K and 
below T=25^30 K, respectively. The dashed lines are T 3 
functions characteristic of c?-wave superconductivity. The 
steep decrease of (1/Ti)host and the upturn of 1/ti is 
also observed for YBa2Cu4_ a; Ni a ;08 with x=0.12 (T c =15 
K) [M. It is theoretically suggested that the spin-orbit 
coupling between an itinerant electron and a Ni local mo- 
ment induces a local superconducting state with a differ- 
ent order parameter (ci^-wave symmetry) around Ni in 
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the da;2_j,2-wave superconducting state [p5{ . The steep 
decrease of (1/Ti)host far below T c may indicate a Ni- 
induced impurity band associated with the different or- 
der parameter, more gapped on the Fermi surface than a 
pure d x 2_ y 2-wave gap. 

Figure 5 shows the temperature dependence of 
(1/Ti)host (a) an d 1/n (b) at Cu(2) in linear scale 
to show up the above T c data. The host (1/Ti)host 
of Cu(2) above T c is nearly independent of Ni doping, 
whereas 1/n increases with decreasing T c . In general, 
1 /ti is an increasing function of the impurity concentra- 
tion ^6|. Thus, the observed increase of 1/n indicates 
the systematic increase of the in-plane Ni(2) concentra- 
tion Xp[ ane . 

In the theoretical model of superconducting pairing 
mediated by antiferromagnetic spin fluctuations, the ex- 
ternal depairing effect on T c is written by T c — T c q — 
AT C with T c0 ~ T a (Q)xa{Q) (T (Q) is the host an- 
tiferromagnetic spin-fluctuation frequency, and Xo(Q) 
is the static staggered spin susceptibility) |3^j3^] and 
with AT C oc Xpi an e fjsj. In the spin- fluctuation the- 
ory §§], one can obtain (I/T^host oc TC/(T + 9) 
(C oc Xo(Q)/r (<9), and 6 oc T (Q) @) in the lead- 
ing order. The actual fit result by this function is the 
dotted curve in Fig. 5(a). Thus, the quantitative tem- 
perature dependence of {1/T{)host tells us the charac- 
teristic spin-fluctuation parameters, Xo{Q) and Tq(Q), 
which may describe the pairing interactions. The nearly 
Ni-independent (1/Ti)host m Fig- 5 indicates that the 
host spin-fluctuation spectrum is nearly invariant un- 
der Ni doping. In Ref. the host {1/Ti)host was 
estimated to be systematically enhanced by Ni doping, 
which was regarded as evidence for softening of the spin- 
fluctuation frequency Fo(<5) and for the central origin 
to reduce T c . However, our analysis indicates that the 
Ni-enhanced Cu(2) nuclear spin-lattice relaxation comes 
from the extra relaxation in the stretched exponential 
part. The central origin to reduce T c is the external pair 
braking effect due to the increase of the in-plane Ni(2) 
concentration x p i ane in AT C . This is consistent with the 
original result for YF^Cu^zNi^Og. However, it is still 
hard to estimate the quantitative value of x p i ane . 

Here, we analyze the temperature dependence of 1/n 
and discuss the pair breaking mechanism of Ni. If an 
impurity spin-spin relaxation process dominates the im- 
purity spin correlation, 1/ti may not change with tem- 
perature |pq| . The widely known expression of 1/Ti or 
1/ti oc {A/h) 2 S{S + l)/uj e (A is the coupling constant of 
the nuclear moment to an electron, S is the electron spin, 
uj e is the exchange frequency of mutual spin interaction) 
is derived from exchange narrowing limit ^Tj.^J. How- 
ever, the estimated 1/n changes with temperature as in 
Figs. 4(b) and 4(d). This indicates a significant contri- 
bution from an impurity spin-lattice relaxation process. 

The temperature dependences of 1/ri far below T c in 
Figs. 4(b) and 4(d) are different between two Ni-doped 



samples. Both the impurity spin-lattice relaxation rate 
Ts l and the impurity spin-spin relaxation rate Tss con- 
tribute to the actual total impurity relaxation rate F^, 
e.g. Tl — Tsl + ^ss- As temperature is decreased, Tsl 
decreases moderately but rapidly below T c 42 and 
then Tss can play an important role in T^. In general, 
Tss is proportional to the number density of impurity 
Ximp but Tsl is independent of Xi mp , and 1/n is pro- 
portional to x 2 mp S(u>N /Tl) (S(lun/^l) is the impurity 
magnetic spectral function at the nuclear resonance fre- 
quency lon /2tt) [^6|,^7|. Thus, the difference in 1/ri may 
result from the different weight of Tss with Xi mp . In 
the below analysis, however, we neglect Tss and then 
the difference in the temperature dependence of 1/ri for 
simplicity. 

For an isolated local moment system in a conventional 
metal, the dynamical spin susceptibility as a function of 
frequency u>/2ir is expressed by x(g, to) — XL /(I— i^/T^), 
where xl oc S(S + 1)/T is the static spin susceptibility 
and T L = aT = 47r( JN F ) 2 k B T /% ( J is the coupling con- 
stant of the localized moments to the band, and Np is the 
density of states at the Fermi level per spin) is the impu- 
rity fluctuation frequency due to Korringa relaxation [ 44 1 . 
Then, one obtains 1/Ti or 1/n oc T l /(uj% + T\) @g7|. 
Since 1/n oc 1/aT at high temperatures (T 3> u>pf/a) 
and 1/n oc aT/ufj at low temperatures (T <C un/cx), 
then 1/n takes a maximum value at T = uj^/a (Tl = 

L0 N ). 

In Figs. 4(b) and 4(d), 1/n of 63 Cu(2) and of 
63 Cu(l) takes a maximum at about 60 K and 10~20 
K, respectively, which can be associated with the maxi- 
mum at T = lvn/ck. Assuming the Korringa relaxation 
T L =4ir(JN F ) 2 k B T/h and putting T L =u N =22 MHz for 
Cu(l) at T=10 K and r L =31.5 MHz for C u(2) at 
T=60 K, one obtains | JN F \=0.b^/Tiuj N /nk B T =0.003 
for Cu(l) and =0.0014 for Cu(2). From a typical 
value of Np=1.5 states/eV-spin direction, |J| is 1.9 
meV for Cu(l) and 0.94 meV for Cu(2). The ac- 
tual maximum of 1/n takes place in the superconduct- 
ing state. Thus, replacing the Korringa term k B T by 
a d-wave gapped term k B T(k B T/ A max ) 2 {2A max = 
Sk B T c HI), one estimates |JTV F |=0.009 and then 
|J|=5.8 meV for Cu(2) (if taking into account the ef- 
fect of the antiferromagnetic spin correlation, one gets a 
smaller |J|). The estimated magnitude of |J| is smaller 
than the in-plane exchange interaction | JAr.;_7Vi| = 11^31 
meV of La 2 Ni0 4 +5 |f||, |J C «-C„|=150 meV of 
YBa 2 Cu 3 06+5 or a typical 4s-3d exchange inter- 

action \J sd \ ~0.1 eV. Since \ttJN f S\ = 0.03 <1, the 
decrease of T c due to the magnetic scattering can be cal- 
culated within the lowest order Born approximation p8[ . 
From AT C = Q.25ir 2 x plane N F J 2 S(S+ l)/k B with 5=1 in 
a d x 2_ y 2-w&ve superconductor |pQ], the sole occupation at 
Cu(2) (x p i ane —x/ '2=0. 05) in the quenched sample yields 
AT C =0.08 K at most. This value is smaller than the ob- 
served AT C ^20 K in the quenched sample by a factor 
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^250, suggesting that magnetic pair breaking is not the 
mechanism of suppression of T c in YBa 2 Cu3_ 2; Ni a ;07_i5 
(similar estimation for Zn doping is seen in Ref. ]5l|). 

Here, one may doubt the assumption of Tl—lon at 
the temperature of maximum ljr\. Instead of 1^=0^, 
one may associate the maximum of 1/ri with a min- 
imum behavior of Tl(^$> ujn) as a crossover from the 
high temperature Tsl to the low temperature T$s- As 
mentioned above, Tss is proportional to x p i ane , so that 
the crossover temperature must increase with increasing 
Xpiane- However, the temperature of maximum 1/ri is 
nearly the same in two samples. Thus, this scenario is 
unlikely. In passing, if the Ni spin freezing temperature 
Tm ~1.5 K (x ^0.09) in Ref. (34) is assigned to the max- 
imum temperature T = uj^/a, one obtains \JNp\=1.9 
and then |J| = 1.3 eV. This magnitude of \ J\ is too large, 
although it can easily reproduce the observed AT C . In 
addition, the Ni spin fluctuation frequency Tl increases 
up to ~20 meV at T c due to the large |J|, so that the 
neutron scattering technique must probe the increase of 
magnetic response around 20 meV due to Ni spin. But, 
the existing data do not support the observation of Ni 
spin fluctuation ]54j ]. Thus, Tm ~ u^/a is also unlikely. 

From AT C = 0.25ir 2 x plane N F J 2 S(S + l)/k B with 
0< x p i ane <0.05 in a d x 2_ y 2-wave superconductor [Q, 
\J\ >68 meV (S=l) or |J| >111 meV (5=1/2) is re- 
quired to account for AT C ^20 K in the quenched sample 
(similar estimation for Zn doping is seen in Refs. f32| , |53| ]). 
Then, \ttNfJS\ >0.3 indicates moderately strong scat- 
tering. The expression of AT C should be extended to 
include both potential scattering and strong scattering. 

In considering the magnetic impurity effect in s-wave 
superconductors, one can neglect nonmagnetic potential 
scattering which should accompany each magnetic impu- 
rity, because of Anderson's theorem, unless the localiza- 
tion effect takes place. However, the situation is com- 
pletely different in the e?-wave state. The potential scat- 
tering also affects d-wave superconductivity as well as the 
magnetic impurity scattering; thus one h as to take into 
account both effects (see also 



. Let us treat 

magnetic impurities as classical spins (S z — ±S) and 
assume that they also have nonmagnetic potential scat- 
terings described by u. Taking into account the impurity 
scatterings within a i-matrix approximation, Ohashi de- 
rived the theoretical expression of T c (x p i ane <Cl) in a 
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y 2-wave superconductor ^ 



\T = "-t" a " e |i 1 ( \ I L_ 

c Ak B N F l 2 V l + ( 7 „- 7m ) 2 l + ( 7 n + 



7m) 2 
(1) 
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Here, "f n and 7m , respectively, describe the depairing ef- 
fects originating from nonmagnetic and magnetic parts 
of the magnetic impurities: 



Putting 7 m = in eq. (1), we obtain AT C in the 
case of nonmagnetic impurity Zn, being in agreement 
with the expression in Refs. [|58[-|6l|]. On the other hand, 
7„ = gives AT C in the case of magnetic impurity with- 
out nonmagnetic potential scattering part. In the limit 
of 7, n — * and 7n — ► 0, eq. (1) leads to the expression 
of AT C in the lowest order Born approximation, being 
in agreement with the expression in Ref. fl50| . In addi- 
tion, 7, n — > oo or 7„ — > co corresponds to the unitarity 
limit |J; AT C = x plane /4k B N F « (1.9 x 10 3 K)x plane . 
If cancellation of "f m — "f n = in the unitarity limit 
takes place for Ni, one obtains AT C = x p i ane /8kBNp, 
which is two times smaller than that due to pure nonmag- 
netic scattering. For YBa 2 Cu 3 _ :E Ni 2 ;07_ ( 5, however, one 
should assume 7„ ^> 7to . In the moderately strong scat- 
tering (a finite u), since AT C oc x p i ane Np at N F <C 1/nu 
(7„ < 1) and AT C /N F at Np > 1/ttu (7„ > 1), 

then AT C takes a maximum of 0.125x p i ane /kBN F at the 
optimal N F = 1/nu ( 7n = 1). If N F of YBa 2 Cu 3 07_ (5 
is optimal, that is u—l/nN F —212 meV (7„ = 1), using 
AT C = (9.5 x 10 2 K)x p iane, we estimate x p / ane =0.01 (as- 
synthesized) and 0.02 (quenched). In the unitarity limit, 
eq. (1) does not explicitly depend on u nor J. Then, 
using AT C = (1.9 x 10 3 K.)x p i ane , we estimate the min- 
imum values of a; p / ane =0.005 (as-synthesized) and 0.01 
(quenched). These values are so reasonable as to sat- 
isfy 0< x p i a ne <0.05. Thus, we suggest that the moder- 
ately strong, nonmagnetic scattering of Ni is a promis- 
ing origin of T c suppression. This is qualitatively con- 
sistent with a report on a large \u\(> | JS\) of Ni in 
Bi 2 Sr 2 CaCu 2 8+d - §§. 

In conclusion, the Cu NQR experiment demonstrated 
that both Cu(l) and Cu(2) nuclear spin-lattice relax- 
ations in YBa 2 Cu3_ a; Ni a ;07_ ( 5 are affected by Ni doping, 
that is, the doped Ni impurities are substituted for both 
Cu(l) and Cu(2) sites. One of the reasons for the small 
decrease of T c by Ni doping in YBa 2 Cu307_^ is partial 
substitution of the Ni impurities for the chain site. In the 
light of the Ni-induced nuclear spin-lattice relaxation, the 
host Cu(2) spin fluctuation spectrum in YBa 2 Cu307_a 
with optimal oxygen content is quite robust for Ni dop- 
ing. 

We would like to thank Professors M. Matsumura, H. 
Yamagata, and K. Mizuno for helpful discussion. This 
work was supported by the New Energy and Industrial 
Technology Development Organization (NEDO) as Col- 
laborative Research and Development of Fundamental 
Technologies for Superconductivity Applications. 



7„ = nN F u, 
7,„ = nN F JS. 



(2) 



[1] G. Xiao, M. Z. Cieplak, J. Q. Xiao, and C. L. Chien, 
Phys. Rev. B 42, 8752 (1990). 



4 



[2] N. Watanabe, N. Koshizuka, N. Seiji, and H. Yamauchi, 

Physica C 234, 361 (1994). 
[3] T. Miyatake, K. Yamaguchi, T. Takata, N. Koshizuka, 

and S. Tanaka, Phys. Rev. B 44, 10139 (1991). 
[4] A. J. Millis, S. Sachdev, and C. M. Varma, Phys. Rev. B 

37, 4975 (1988). 
[5] R. Balian, and N. R. Werthamer, Phys. Rev. 131, 1553 

(1963). 

[6] J. T. Markert, B. D. Dunlap, and M. B. Maple, MRS 

Bulletin, 14, 37 (1989), and the references therein. 
[7] E. Takayama-Muromachi, Y. Uchida, and K. Kato, Jpn. 

J. Appl. Phys. 26, L2087 (1987). 
[8] J. M. Tarascon, L. H. Greene, P. Barboux, W. R. McK- 
innon, G. W. Hull, T. P.Orlando, K. A. Delin, S. Foner, 
and E. J. McNiff, Jr., Phys. Rev. B 36, 8393 (1987). 
[9] Y. Tokunaga, K. Ishida, Y. Kitaoka, and K. Asayama, 
Solid State Commun. 103, 43 (1997). 

[10] T. Kluge, Y. Koike, A. Fujiwara, M. Kato, T. Noji, and 
Y. Saito, Phys. Rev. B 52, 727 (1995). 

[11] Y. Itoh, T. Machi, N. Watanabe, S. Adachi, and N. 
Koshizuka, J. Phys. Soc. Jpn. 70, 1881 (2001). 

[12] Y. Itoh, T. Machi, N. Watanabe, and N. Koshizuka, 
Physica C 357-360, 69 (2001). 

[13] K. Ishida, Y. Kitaoka, N. Ogata, T. Kamino, K. 
Asayama, J. R. Cooper, and N. Athanassopoulou, J. 
Phys. Soc. Jpn. 62, 2803 (1993). 

[14] R. S. Howland, T. H. Geballe, S. S. Laderman, A. 
Fischer-Colbrie, M. Scott, J. M. Tarascon, and P. Bar- 
boux, Phys. Rev. B 39, 9017 (1989). 

[15] C. Y. Yang, A. R. Moodenbaugh, Y. L. Wang, Youwen 
Xu, S. M. Heald, D. O. Welch, M. Suenaga, D. A. Fischer, 
and J. E. Penner-Hahn, Phys. Rev. B 42, 2231 (1990). 

[16] F. Bridges, J. B. Boyce, T. Claeson, T. H. Geballe, and 
J. M. Tarascon, Phys. Rev. 42, 2137 (1990). 

[17] K. Morita, H. Yamagata, and M. Matsumura, in Meet- 
ing Abstracts of the Physical Society Japan, Higashi- 
Hiroshima, 1999, Vol. 54, p. 559 [in japanese]. 

[18] S. Adachi, Y. Itoh, T. Machi, E. Kandyel, S. Tajima, and 
N. Koshizuka, Phys. Rev. B 61, 4314 (2000). 

[19] C. C. Homes, D. A. Bonn, R. Liang, W. N. Hardy, D. N. 
Basov, T. Timusk, and B. P. Clayman, Phys. Rev. B 60, 
9782 (1999). 

[20] T. Wada, N. Suzuki, T. Maeda, A. Maeda, S. Uchida, K. 
Uchinokura, and S. Tanaka, Appl. Phys. Lett. 52, 1989 
(1988). 

[21] S. I. Yoo, N. Sakai, H. Takaichi, T. Higuchi, and M. Mu- 
rakami, Appl. Phys. Lett. 65, 633 (1994). 

[22] S. Shibata, A. K. Pradhan, and N. Koshizuka, Jpn. J. 
Appl. Phys. 38, L1169 (1999). 

[23] A. Navrotsky and O. J. Kleppa, J. Inorg. Nucl. Chem. 
29, 2701 (1967). 

[24] Y. Itoh, T. Machi, N. Watanabe, and N. Koshizuka, J. 
Phys. Soc. Jpn. 68, 2914 (1999). 

[25] Y. Itoh, T. Machi, N. Watanabe, and N. Koshizuka, J. 
Phys. Soc. Jpn. 70, 644 (2001). 

[26] M. R. McHenry, B. G. Silbernagel, and J. H. Wernick, 
Phys. Rev. Lett. 27, 426 (1971); Phys. Rev. B 5, 2958 
(1972). 

[27] T. Moriya, Prog. Theor. Phys. 16, 23 (1956); 16, 641 

(1956); 28, 371 (1962). 
[28] C. H. Pennington, D. J. Durand, D. B. Zax, C. P. 



Slichter, J. P. Rice, and D. M. Ginsberg, Phys. Rev. B 

37, 7944 (1988). 
[29] T. Shimizu, H. Yasuoka, T. Imai, T. Tsuda, T. Taka- 

batake, Y. Nakazawa, and M. Ishikawa, J. Phys. Soc. 

Jpn. 57, 2494 (1988). 
[30] J. Bobroff, H. Alloul, Y. Yoshinari, A. Keren, P. Mendels, 

N. Blanchard, G. Collin, and J.-F. Marucco, Phys. Rev. 

Lett. 79, 2117 (1997). 
[31] G. R. Khutsishvili, Soviet Phys. Uspekhi, 8, 743 (1966). 
[32] G. V. M. Williams, and S. Kramer, Phys. Rev. B 64, 

104506 (2001). 

[33] E. Kukovitsky, H. Liitgemeier, and G. Teitel'baum, Phys- 
ica C 252, 160 (1995). 

[34] Y. Tokunaga, K. Ishida, Y. Kitaoka, and K. Asayama, 
Czech. J. Phys. 46, 1139 (1996) S2. 

[35] A. V. Balatsky, Phys. Rev. Lett. 80, 1972 (1998). 

[36] T. Moriya, and K. Ueda, J. Phys. Soc. Jpn. 63, 1871 
(1994). 

[37] P. Monthoux, and D. Pines, Phys. Rev. B 49, 4261 
(1994). 

[38] A. A. Abrikosov, and L. P. Gorkov, Sov. Phys. JETP 12, 
1243 (1961). 

[39] T. Moriya, Y. Takahashi, and K. Ueda, J. Phys. Soc. 

Jpn. 59, 2905 (1990). 
[40] Y. Itoh, Physica C 263, 378 (1996). 
[41] P. W. Anderson, and P. R. Weiss, Rev. Mod. Phys. 25, 

269 (1953). 

[42] M. Mali, D. Brinkmann, L. Pauli, J. Roos, H. Zimmer- 
mann, and J. Hullinger, Phys. Lett. A 124, 112 (1987). 

[43] T. Imai, T. Shimizu, H. Yasuoka, Y. Ueda, and K. Ko- 
suge, J. Phys. Soc. Jpn. 57, 2280 (1988). 

[44] A. C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, 1997). 

[45] N. Bloembergen, E. M. Purcell, and R. V. Pound, Phys. 
Rev. 73, 679 (1948). 

[46] T. Hasegawa, H. Ikuta, and K. Kitazawa, in Physical 
Properties of High Temperature Superconductors 
III, edited by D. M. Ginsberg (World Scientific, Sin- 
gapore, 1992). 

[47] K. Yamada, M. Arai, Y. Endoh, S. Hosoya, K. Nakajima, 
T. Perring, and A. Taylor, J. Phys. Soc. Jpn. 60, 1197 
(1991). 

[48] K. Nakajima, Y. Endoh, K. Yamada, S. Hosoya, N. 
Kaneko, and J. Wada, J. Phys. Soc. Jpn. 64, 716 (1995). 

[49] J. Rossat-Mignod, L. P. Regnault, C. Vettier, P. Bourges, 
P. Burlet, J. Bossy, J. Y. Henry, and G. Lapertot, Physica 
C 185-189, 86 (1991). 

[50] S. Maekawa, Y. Isawa, and H. Ebisawa, Jpn. J. Appl. 
Phys. 26, L771 (1987). One should note that this expres- 
sion of AT C in the lowest order Born approximation for 
a d-wave superconductor is smaller than the Abrikosov- 
Gorkov formula AT C = 0.5n 2 x imp N F J 2 S(S + l)/k B for 
an s-wave superconductor by a numerical factor 2. See 
also L. A. Openov, Phys. Rev. B 58, 9468 (1998). 

[51] R. E. Walstedt, R. F. Bell, L. F. Schneemeyer, J. V. 
Waszczak, W. W. Warren, Jr., R. Dupree, and A. 
Gencten, Phys. Rev. B 48, 10646 (1993). 

[52] A. V. Mahajan, H. Alloul, G. Collin, and J. F. Marucco, 
Phys. Rev. Lett. 72, 3100 (1994). 

[53] S. Zagoulaev, P. Monod, J. Jegoudez, Physica C 259, 
271 (1996). 



5 



[54] Y. Sidis, P. Bourges, H. F. Fong, B. Keimer, L. P. 

Regnault, J. Bossy, A. Ivanov, B. Hennion, P. Gautier- 

Picard, G. Collin, D. L. Millius, and I. A. Aksay, Phys. 

Rev. Lett. 84, 5900 (2000). 
[55] M. I. Salkola, A. V. Balatsky, and J. R. Schrieffer, Phys. 

Rev. B 55, 12648 (1997). 
[56] R. Kilian, S. Krivenko, G. Khaliullin, and P. Fulde, Phys. 

Rev. B 59, 14432 (1999). 
[57] Y. Ohashi, unpublished works. 

[58] L. S. Borkowski, and P. J. Hirschfeld, Phys. Rev. B 49, 
15404 (1994). 

[59] R. Fehrenbacher, and M. R. Norman, Phys. Rev. B 50, 
3495 (1994). 

[60] P. Arberg, and J. P. Carbotte, Phys. Rev. B 50, 3250 
(1994). 

[61] Y. Kitaoka, K. Ishida, and K. Asayama, J. Phys. Soc. 

Jpn. 63, 2052 (1994). 
[62] T. Hotta, J. Phys. Soc. Jpn. 62, 274 (1993). 
[63] E. W. Hudson, K. M. Lang, V. Madhavan, S. H. Pan, 

H. Eisaki, S. Uchida, and J. C. Davis, Nature 411, 920 

(2001). 




YBa 2 (Cu,_ z M z ) 3 7 
(a) 



Bi 2 Sr 2 Ca(Cu,_ z M z ) 2 8 - 
YBa 2 (Cu,_ z M z ) 4 8 



\ V==— ' La i.85 Sr o 15CU1.MO4 



YBa 2 (Cu, _MJ 2 7 

N^as-synthesized" M 




0.05 0.10 
M(=Ni, Zn) concentration z 



0.00 0.01 0.02 0.03 
Ni concentration z 



FIG. 1. (a)Impurity doping dependence of T c as func- 
tions of the concentration z for various high-T c supercon- 
ductors with M=Ni or Zn. The solid (open) symbols 
are the Ni(Zn) doping dependence of T c . The data are 
adopted from [9] for YBa 2 (Cui_ z M z ) 3 7 (Y123, upward tri- 
angles), [10] for Bi 2 Sr 2 Ca(Cui_ z M z ) 2 8 (Bi2212, squares), 
[11,12] for YBa 2 (Cui_ z M z ) 4 8 (Y124, circles), and [1] for 
Lai.85Sro.i5Cui_ z M z 04 (LS214, downward triangles). The 
dashed and the solid lines are fit by T c =T c o-rriMZ (niM=Ni,Zn 
is the fitting parameter) for the respective materials. The esti- 
mated ratio rriNi/mzn is about 0.26 for Y123, 0.46 for Bi2212, 
0.80 for Y124, and 0.62 for LS214. The decrease of T c by Ni 
doping for Y123 is smallest among these materials. (b)Ni 
doping dependence of T c for "as- synthesized" or "quenched" 
Y123, adopted from [18]. The solid curves are guide for the 
eyes. Note that the Ni concentration x in the text is denned 
bya;=3ainYBa 2 Cu3- :c M x 07_ ( 5 (7-5=6.92-6.95 [18]). The de- 
crease of T c by Ni doping is larger in "quenched" Y123 than 
in "as-synthesized" one. 
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FIG. 2. Zero-field frequency spectra of the chain-site and 
the plane-site 63 ' 65 Cu NQR for YBa 2 Cu 3 -^Ni ;c 07-i at T=4.2 
K. The line shapes are scaled after T 2 corrections have been 
made. 
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FIG. 3. The Ni-doping effect on the 63 Cu nuclear 
spin-echo recovery curves p(t) = 1 — M(t)/M(po) at T=4.2 K. 
Inset figures show the recovery curves for the pure (Ni-free) 
sample at T=4.2 K. The solid curves are the least-squares 
fitting results using the theoretical function including a 
stretched exponential function (see the text). 
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FIG. 4. Log-log plots of 63 (1/Ti)host (a) [(c)] and 
63 (l/n) (b) [(d)] at 63 Cu(2) [ 63 Cu(l)] as functions of tem- 
perature for the Ni-free and for the Ni-doped samples. The 
dashed lines are T 3 functions. The solid curves are guide for 
the eyes. 
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FIG. 5. The Ni-doping effect on the temperature depen- 
dence of 63 (1/Ti) H ost (a) and 63 (1/t!) (b) at 63 Cu(2) in 
linear scale. The thick bars indicate the respective T c 's. The 
dotted curve is the least-squares fitting result by a function 
of TC/(T + 9) ((7=1.7 ms- 1 and 9=128 K) above T c . 
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